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Abstract. We present an improved algorithm for the computation of Zariski cham- 
bers on algebraic surfaces. The new algorithm significantly outperforms the so far 
available method and allows therefore to treat surfaces of high Picard number, where 
huge chamber numbers occur. As an application, we efficiently compute the number 
of chambers supported by the lines on the Segre-Schur quartic. 



Introduction 

Zariski chambers are natural pieces into which the big cone of an algebraic surface 
decomposes. Their properties were first studied in It is an intriguing problem, 
raised in [3], to determine explicitly how many Zariski chambers a given surface has. 
In other words, we ask on a smooth projective surface X for the quantity 

z(X) =def # {Zariski chambers on X} £ N U {oo} . 

Roughly speaking, the number z(X) measures how complicated the surface X is 
from the point of view of linear series. Specifically, it answers the following natural 
questions (see [3]): 

• How many different stable base loci can occur in big linear series on XI 

• How many essentially different Zariski decompositions can big divisors on X 
have? (Here we consider Zariski decompositions to be essentially different, if 
their negative parts have different support.) 

• How many "pieces" does the (piecewise polynomial) volume function on the 
big cone of X have? 

To get a more detailed picture of the geometry of X, it is also natural to consider 
for integers d ^ 1 the numbers 

, . J Zariski chambers on X that are 

dcf \ supported by curves of degree ^ d 

where the degree of curves is taken with respect to a fixed ample line bundle on X. 
One has of course z{X) = sup d z^{X). While there are surfaces X where z(X) is 
infinite, the numbers Zd{X) are always finite, because on any smooth surface there 
are only finitely many irreducible negative curves of fixed degree. 
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In [3] an algorithm was presented that computes Zariski chambers from the 
intersection matrix of a set of negative curves, and the algorithm was applied to Del 
Pezzo surfaces. While this method is very efficient in these cases, further experience 
has shown that there do exist surfaces where the algorithm takes an inordinate 
amount of time - to the point of becoming impractical in such situations. This 
is for instance the case when the method is being applied to the 64 lines on the 
Segre-Schur quartic (see Sect. [2]). At first this phenomenon is somewhat surprising, 
as there are surfaces with many more curves to be considered (like the Del Pezzo 
surface X$ with its 240 negative curves) where an application of the algorithm poses 
no practical problems at all. It seems that it is not so much the number of negative 
curves that matters most, but the Picard number of the surface and the number 
of chambers that are found. This is in accordance with the observation that the 
essential work done by the algorithm (and in fact its essential bottle-neck) is the 
computation of an enormous number of determinants - and the dimension of the 
matrices in question is bounded in terms of the Picard number of the surface. (The 
dimension of the determinants to be computed is bounded by p(X) — 1.) 

In the present paper we attack this problem by providing a significantly improved 
algorithm that is suited also for surfaces with higher Picard number. Our focus here 
is on the efficient calculation of the relevant determinants. As is well-known, usual 
fraction-free algorithms for the computation of the determinant of an n x n integer 
matrix, like the fraction-free Bareiss algorithm, have complexity 0(n 3 ). And, to our 
knowledge, even the currently best fraction-free algorithms for determinants over 
integral domains need (9( n 2 - 697263 ) r i n g operations (see [SJ). Our main point is that 
within our new algorithm each of the necessary determinant computations has a 
complexity of only 0(n 2 ) - this is achieved by reusing information from previous 
computations (see the details in Sect. [1]). 

As an application, we determine the number of chambers that are supported 
by the lines on the Segre-Schur quartic. This remarkable surface provides an ideal 
application for the new algorithm: It turns out that it has an enormous number of 
Zariski chambers supported by lines, and it seems that the surface lies at the edge 
of what can be practically computed with current methods. 

We show: 

Theorem. Let X C P 3 be the Segre-Schur quartic, i.e., the surface given in homo- 
geneous coordinates {xq : x% : X2 ■ £3) by the equation 

x (xl - x\) - x 2 {x\ - x\) = . 

(i) We have 

z(X) = 00 

and 

Zi(X) = 8 260 383 569. 

(ii) The maximal number of lines that occur in the support of a Zariski chamber 
is 19 (which is the maximal possible value, as the Picard number of X is 20). 

Note that the number z±(X) is bigger by a factor of about 10 5 than the number 
1 501 681 of chambers on the Del Pezzo surface X$ (blow-up of the plane in 8 points) - 
even though only 64 curves are used to build chambers, as opposed to 240 curves 
on X$. On the other hand, the Picard number of the Segre quartic is about twice 
that of X$. One is lead to wonder in general how the number of Zariski chambers is 
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related to the Picard number - apart from the fact that with higher Picard number, 
chambers of bigger support size become theoretically possible. 

Note that since the negative curves on X of higher degrees are not known, the 
numbers Zd(X) are at present inaccessible for d ^ 2. It would already be very 
interesting to know at what rate they grow for d — > oo. 

Concerning the organization of this paper, we start in Sect. 1 by presenting the 
improved algorithm. In Sect. 2 we give the proof of the theorem on the Segre-Schur 
quartic stated above. Finally, Sect. 3 compares the new algorithm with the original 
one, providing sample run-times for the case of Del Pezzo surfaces, the Segre-Schur 
quartic, and for matrices related to Fermat surfaces. 



Zariski chambers were first studied in [3J, and we refer to that paper for a detailed 
introduction. For a very brief account, consider a smooth projective surface X over 
the complex numbers. To any big and nef ]R-divisor P on X, one associates the 
Zariski chamber £p, which by definition consists of all divisor classes [D] in the 
big cone Big(X) such that the irreducible curves in the negative part of the Zariski 
decomposition of D are precisely the curves C with P ■ C = 0. It is the main result 
of [3J that the sets £p yield a locally finite decomposition of Big(X) into locally 
polyhedral subcones such that 

• on each subcone the volume function is given by a single polynomial of degree 



• in the interior of each of the subcones the stable base loci are constant. 

For the purposes of the present paper it is a crucial fact that the number of Zariski 
chambers can be computed from the intersection matrix of the negative curves on X, 
because Zariski chambers correspond to negative definite reduced divisors - one has 



Proposition 1.1. The set of Zariski chambers on a smooth projective surface X 
that are different from the nef chamber is in bijective correspondence with the set of 
reduced divisors on X whose intersection matrix is negative definite. 

In the statement, the nef chamber is the chamber £p associated with an ample 
divisor H. Its closure is the nef cone, and its interior is the ample cone. Suppose 
next that X carries only finitely many negative curves (i.e., irreducible curves C C X 
with negative self-intersection C 2 ). Then one has, as an immediate consequence of 
Prop. O (see [3, Prop. 1.5]): 

Proposition 1.2. (i) The number z(X) of Zariski chambers on X is given by 



(ii) More generally, let C\, . . . ,C r be distinct negative curves on X, and let S 
be their intersection matrix. Then the number of Zariski chambers that are 
supported by a non-empty subset of {C\, . . . , C r } equals the number of negative 
definite principal submatrices of the matrix S. 



1. Efficient computation of Zariski chambers 



two, and 



by [31 Prop. 1.1]: 




negative definite principal submatrices of the 
intersection matrix of the negative curves on X 
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Here a principal submatrix of a given (n x ra)-matrix is as usual a submatrix that 
arises by deleting k corresponding rows and columns of the matrix, where ^ k < n. 

The subsequent algorithm computes the number of positive definite principal 
submatrices of a given symmetric matrix. In view of Proposition 11.21 this enables 
us to determine the number of Zariski chambers (by considering the negative of the 
intersection matrix). 

Algorithm 1.3. The algorithm takes as input an integer n 1 and a symmetric 
(n x n)-matrix A over Z. It outputs all subsets S C {1, . . . , n} having the property 
that the corresponding principal submatrix As is positive definite. 

We adopt a backtracking strategy as in [3], but instead of testing for positive 
definiteness by a computation of the determinant det(-As) from scratch, the algo- 
rithm makes use of three procedures Grow, Shrink, and IsPosDef to be discussed 
below. 

Algorithm PosDef 

Input: integer n 1, symmetric matrix A £ Z nxn 
Output: all positive definite principal submatrices of A 

k <- 1; S <- 0; B <- (); T <- () 
Grow(S, k) 
while S 7^ do 

Assert {{B = T ■ As) and {B is in Bareiss form)) 

if IsPosDef (5) then 
output S 

else 

Shrink (S) 
end if 

if k < n then 

k <— k + 1; Grow(S,k) 

else 

if S 7^ and A; = max 5 then 

Shrink(S) 
end if 

if S / then 

maxS; Shrink(S); k •<— A; + 1; Grow(S, k) 

end if 
end if 
end while 

We now explain the procedures Grow, Shrink, and IsPosDef. They work on 
matrix variables B and T that are global variables of PosDef. At every stage of the 
algorithm, B holds the Bareiss trigonalization of As- The auxiliary matrix variable 
T is the essential tool that makes is possible to do the trigonalization incrementally. 
It holds the triangular matrix that one obtains from the unity matrix, when the 
same transformations are applied that have been applied As to obtain B. 

The procedure Shrink(S) removes the maximal element from S and updates B 
and T. The latter is achieved by simply discarding the last row and column of both 
B and T. 

The procedure IsPosDef (S) determines whether As is positive definite. This can 
be done as follows: As the matrix A5"\ max 5 is known to be positive definite, As is 
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positive definite if and only if its determinant is positive. And the latter can be read 
off the sign of the lower right entry of B, since in a Bareiss trigonalization this entry 
is always the determinant of the original matrix As- So we have: 

IsPosDef(S) = (S[max5,max5] > 0) . 

The procedure Grow(S, k) includes a new element k > max S to the index set S 
and updates B and T accordingly. It consists of the following steps: 

(Gl) S^Su{k} 

(G2) Attach the last row and the last column of As to the bottom and right of B. 
Attach the last row and column of the unit matrix to T. 

(G3) Replace the last column b of B by T ■ b . 

(G4) Clear the last row of B by means of the following procedure: 

s^\S\ 

for i from 1 to s — 1 do 
if i = 1 then 
d<-l 

else 

d<- B[i- l,i- 1] 
end if 

for j from i + 1 to s do 

B[s,j] <- (B[s,j]-B[i,i] 
end for 

for j from 1 to s do 

T[s,j] ir- (T[s,j]-B[i,i] 
end for 
B[s,i] «- 
end for 

In the two loops of (G4), the division by d is the Bareiss divison, made possible by 
Sylvester's identity, which keeps the size of the matrix entries from exploding (cf. [lj). 
Note that Grow does 0(s 2 ) loop iterations and hence has complexity 0(n 2 ), while 
Shrink and IsPosDef need only constant time. 

2. Lines and chambers on the Segre-Schur quartic 

In this section we consider the Segre-Schur quartic (see [TU], [12], and also [2J 
Sect. 2.1]). By Segre's theorem [11], this remarkable surface carries the maximal 
number of lines that is possible for a smooth quartic in P 3 . In order to find the 
chambers supported on lines on this particular surface we need to determine the 
intersection matrix of all lines. We approach this task from a more general setting. 
Caporaso, Harris, Mazur [7J and Boissiere, Sarti [5] consider a class of surfaces in P 3 
known to contain many lines, which was first studied by Segre [12] . Namely, these 
surfaces are given by an equation 



- B[i,j] ■ B[s,i]) div d 
-T[i,j}-B[s,i}) divd 



ip(x ,xi) = i/j(x 2 ,x 3 ), 



(2.0.1) 
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with homogeneous polynomials p and ip of common degree d. Clearly the surface 
we are interested in, the Segre-Schur quartic, is of this type with 



cp(x,y) = tp{x,y) = x(x 3 - y 3 ). 

For any surface S given by an equation (|2.0.ip of degree d consider the sets of zeros 
V{tp) and V(ip) in P 1 . Denote by T the group of automorphisms of P mapping 
V{tp) onto Vty). 

Proposition 2.1 ([7j Lemma 5.1, [5j Proposition 4.1). The number of lines 
on S equals d 2 + d ■ \T\ . 

In order to establish terminology for our further investigation, we briefly go 
through the steps of the proof given in [7] . 

• Consider the sets of points Pi , . . . , P^ and P[ , . . . , P' d of S lying on the lines 
L\ = V{x2,xj) and L2 = V{xq,x\) respectively. For every i,j £ {1, ...,d} 
the line L^j joining Pi and P- is contained in S. The d 2 lines Ljj are called 
lines of the first type. 

• Any further line L on S is called a line of the second type. Any such line is 
disjoint from L\ and L2, guaranteeing that it is given by equations of the form 

X2 = oxq + bx\ 
X3 = cxo + dx\ 

with ad — be 7^ 0. The invertible matrix ^ a induces an automorphism of 
P 1 mapping V(tp) onto V(^). 

• Conversely, every automorphism in T given by an invertible matrix 



a b 

y c d 

corresponds to d distinct lines of the second type on S given by equations 



X2 = Xr/ k (axo + bx\) 

x 3 = \ri k {cx Q + dxi), (2.1.1) 
for some A G C and 1 ^ k ^ d where rj denotes a primitive d-th root of unity. 

With this knowledge we turn to explicitly determining the lines lying on the Segre- 
Schur quartic. The 16 lines of the first type are just the lines joining the zeros of 
the polynomial tp(x,y) = x(x 3 — y 3 ) on L± and L2 each considered as a P 1 . Setting 

2Tvi 

£ = e 3 , these points are 

Pi = (0 : 1 : : 0) P 2 = (1 : 1 : : 0) 
P 3 = (£ : 1 : : 0) P 4 = (£ 2 : 1 : : 0) 



and 



P[ = (0 : : : 1) P' 2 = (0 : : 1 : 1) 
P^ = (0 : : e : 1) = (0 : : £ 2 : 1). 
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The lines L 



P{Pj joining them are given as 



Li j : 



-1 dj 






-1 aj 



( x \ 

\ ^ J 



where Oj denotes the i-th entry of the tuple (0, 1,£,£ 2 ). 

For the lines of the second type we note that in the case of the Segre-Schur 
quartic T is the tetrahedral group T, which is isomorphic to the product Di x C3 of 



the dihedral group D2 = (Z/2Z) and a cyclic group C3 = 
the group D2 operating on the points (0 : 1),(1 : 1), (1 
elements 



Z/3Z. More concretely, 
: £)>(! : £ 2 ) consists of 



1 

1 

-1 e 

2£ 1 



-1 £ 

2e 1 

-1 1 
2 1 



and the cyclic group C3 is generated by the element 

Z = 

Hence the elements of T are the automorphisms Z k Tj for ^ k ^ 2 and 



£ 
1 



1 ^ j ^ 4. For an element 7 6 T represented by a matrix 
the corresponding numbers X m such that the lines 



a b 
c d 



we need to find 



L 



X m a X m b 
X m c X m d 



-1 







Xx 
V ^3 / 



from equation (|2.1.ip lie on S. If (p(po,Px) = for a p G P , then 
(^(Aapo + Afcpi, Acp + Adpi) = AV( a Po + bpx,cp + dpi 



0. 



The line L 1) \ thus intersects 5 in the four points (po : px ■ X(apo+bpx) : A(cpo+dpi)), 
corresponding to the zeros (po : Pi) of <p. Now, let (go : <?i) G P 1 be any point on 
which tp does not vanish. Then the desired values for A are given as solutions of the 
equation 

A 4 _ V?(gO,gl) 

ip(aq + bqx,cq + dqx)' 

By way of example, let us carry out the computation for the automorphism 
Z l T2 G T. It is given by the matrix 

-e e 
2e 1 

Choosing (qo : qx) = (1 : 0), we get 

vu, 0) = 1, ^(-e, 2c 2 ) = -e(-e 3 - se 6 ) = 9$. 
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Consequently, the corresponding values for A m are 

Xm-^—, 

for m = 0, . . . , 3. 

In an analogous manner we find the factors X m for the remaining automorphisms 
in T as shown in Table [TJ We set n = v3/3 and write i = \J — 1 for the imaginary 
unit. 



automorphism matrix 



Z°Ti 



Z°Tr 



Z°T 4 



1 
1 



Z-7'j ( C 1 1 



■i e 



1 i 

2 1 



\ 2e i 1 / /; 

( ~t f ) * m ^ 2 

^2 ( ~f 2 \ ) i m vH 

71)-/- / _ 1 \ ;,„, 



^ ( ~l \ ) v-^e 



i m rj 



Table 1: Factors A m for the automorphisms in T 

The lines L Z k T . \ m for A; = 0, . . . , 2, j = 1, . . . , 4, m = 0, . . . , 3 are thus exactly 
the 48 lines of the second type on S. 

The configuration of the lines on the Segre-Schur quartic is given by the following 



Proposition 2.2. Let 6 : R — > {0, 1} denote the indicator function of the set {0} 

mapping any non-zero number x to and the number to 1. VFii/i further notation 
as above, we have: 

• Every line L on S has self-intersection L 2 = —2. 
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The intersection number of distinct lines of the first type is 



L 



S(i-i') + 5(j-j'). 



The intersection number of a line of the first type with a line of the second 
type is 



L s ,fL Z k T • x „ 



5(t + 2k-j-l mod 3) 
5(s — 1 — j mod 3) 
S(s + k — t mod 3) 
[5(s-l)-S(j-l)-S(t-l) 



ifj^l,s^l,t^l, 
ifj^l,t^l,s = l, 

else 



The intersection number of distinct lines of the second type is 

' s^m'-m _ 3 )£2fc 

I jm'—mf^ 3i m ' ~ m ^2fc' 

-3,{i 2m ' v 2 e k ')) , ifj 



L Z k T ] ,\ m - L Z k 'T,,,\ 



5(m 



m 



, if 3= f 



1. 



Proof. The statement about self-intersection of lines on S follows immediately by 
adjunction and the fact that S has trivial canonical class. To prove the statement 
about the intersection of distinct lines of the first type, all we need to show is that 
two such lines A, B cannot intersect in a point outside the (disjoint) lines L\ and L2 
(see sketch of proof of Proposition 12, ip . But if this were the case, since L\ and L2 
both intersect A and B, they would lie in the plane spanned by A and B and would 
therefore intersect each other. 

The remaining intersection numbers are calculated in the following way: Two 
lines in P 3 given by equations 

a x + a\xi + a 2 x 2 + 03^3 = 
b x + bixi + b 2 x 2 + b 3 x 3 = 



and 



CqXq + C\X\ + C 2 X 2 + C3X3 = 

d x + d\xi + d 2 x 2 + d 3 x 3 = 



intersect each other if and only if the determinant of the matrix 

\ 



/ a ai a 2 a 3 

bo bi b 2 b 3 

Co ci c 2 c 3 

\ d di d 2 d 3 



J 



vanishes. Note that for s, t 7^ 1 the line L s j of the first type is given by 



-1 £*" 2 






-1 e~ 2 



( X Q \ 

x\ 

X2 
\ ^ J 



0, 
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Figure 1: Intersection matrix of the 64 lines on the Segre-Schur quartic. The numerical entries of 
—2, and 1 are replaced by black, grey and white boxes respectively. 



and if j ^ 1, then the line L Z k T . A is given by 



L 



Z Tj ,A m 



-i m r] 



-1 



V ^3 J 



The computation of the determinants in question then yields the asserted formulas. 
Let us by way of illustration calculate the intersection number of the line Li t of the 
first type, where t / 1, and a line L Z k T . Xm of second type, with j ^ 1. We have 



det 



-1 







\ 

o -l e~ 2 

i m r] i 171 ^- 1 -1 

\ 2i m r)£ 2k+1 -i i m v £ 2k -I J 



-2k+t-2\ 



Consequently, the lines intersect if and only if 



which is equivalent to 



t + 2k-j-l = mod 3. 



□ 



Combining all the intersection numbers yields the intersection matrix of the 
Segre-Schur quartic displayed in Figure [TJ The order of rows and columns was 
chosen as follows: The first 16 rows correspond to the lines of the first type L s t 
with s = 1, . . . , 4, t = 1, . . . , 4. The remaining 48 rows correspond to the lines of the 
second type L zkT . iXm with k = 0, . . . , 2, m = 0, . . . , 3, j = 1, . . . , 4. 
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Proof of the theorem stated in the introduction. The intersection matrix of the 64 
lines is of rank 20, hence the classes of the lines generate the Neron-Severi group 
NS(X) over Q, and the K3 surface is singular (in the sense that NS(JT) is of maximal 
possible rank). Therefore, by a result of Shioda and Inose [13] . its automorphism 
group is infinite. The latter fact in turn implies that there are infinitely many 
(— 2)-curves on X (cf. [9, Remark 7.2]), hence z(X) = oo. 

An application of Algorithm 11.31 to the negative of the intersection matrix of the 
64 lines shows that it has exactly 8 260 383 568 positive definite principal submatri- 
ces. By taking the additional nef chamber into account we arrive at the claimed 
8 260 383 569 chambers. As p(X) = 20, there cannot be Zariski chambers supported 
by more than 19 curves. The results of the algorithm show that chambers with 19 
supporting curves actually occur. (In fact there are exactly 1728 such chambers, see 
the subsequent remark \2. 31 ) □ 

Remark 2.3. There is more numerical data that is of interest: How many Zariski 
chambers are there with support of given cardinality? Let us denote for £ ^ 1 by 
z j (X) the number of Zariski chambers that are supported by I curves of degree d. 
Our computations yield for the Segre-Schur quartic the values in the following table. 



I zf\x) 



1 


64 


2 


2016 


3 


41376 


4 


605856 


5 


6343776 


6 


45613512 


7 


217025520 


8 


674047818 


9 


1376161536 


10 


1900843848 


11 


1832006112 


12 


1264421472 


13 


635795760 


14 


233619648 


15 


61499712 


16 


11037702 


17 


1246368 


18 


69744 


19 


1728 



The numbers show that the surface clearly favours chambers of "medium size" . 
Remark 2.4. The lines 

Ll t l, • • • , il,4, £2,1, -^2,2) -^2,3; £3,1) -^3,2) -^3,3 

of the first type together with the lines 



T2A0 ' 
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of the second type generate a sublattice A C NS(X) of rank 20 and discriminant —48. 
So they yield a basis of NS(X) (g) Q. It is shown in [6j Appendix B] that the 
discriminant of the lattice NS(AT) is —48 as well, hence the mentioned lines in fact 
generate NS(X). 

Remark 2.5. It is interesting to note that the 16 lines of the first type that one 
has on any quartic surface of type 

(p(x ,xi) = Tp(x 2 ,x 3 ) 

give rise to 6521 Zariski chambers. As these 16 lines generate a sublattice of rank 
10, one should compare the number 6521 with the number z(Xg) = 1501681 for 
the Del Pezzo surface X% of Picard number 9. It would be interesting to know on a 
more conceptual basis what properties of the lattice are crucial for obtaining a large 
number of chambers. 



3. Comparison of efficiency 

In order to illustrate the practicability of the new algorithm, we will now compare 
it with the algorithm from [3] by providing sample run-timesQ 

We consider first the Del Pezzo surfaces X\, . . . , X$ as in [3J. The following table 
lists their chamber numbers z(X r ), and it shows the run-times (all in milliseconds) 
for the algorithm from [3], called A\ in the table, and for the new algorithm A 2 . The 
value n is the dimension of the intersection matrhJ! (i.e., the number of negative 
curves on X r ). 



6 



7 



8 



n 
z(X r 



3 
5 



6 

18 



10 
76 



16 27 
393 2 764 



56 
33 645 



240 
1501681 



Ax 
A 2 



0.55 
0.8 



0.56 0.7 
0.85 0.82 



0.85 
0.91 



2.3 
2.3 



33.5 
26.8 



1.48 x 10 3 2.88 x 10 4 
1.09 x 10 3 1.92 x 10 4 



Factor 0.69 0.66 0.85 0.93 1 1.25 



1.36 



1.5 



For r ^ 4 the original algorithm A± is actually faster, presumably due to the overhead 
caused by the more sophisticated strategy of A 2 . Starting with r = 6, however, A 2 
gets more and more superior. This pattern will become even more visible when we 
now consider the lines on the Segre quartic. Specifically, we provide run-times for 
A\ and A 2 , when applied to the principal submatrices consisting of the first n rows 
and columns of the intersection matrix from Sect. [2 



16 



24 



32 



40 



56 



64 



Ai 
.4, 



0.9 
0.96 



65,55 
31.54 



4.16 x 10 J 
1.7 x 10 3 



9.46 x 10 4 
3.61 x 10 4 



1.43 x 10 b 
5.4 x 10 5 



1.57 x 10 7 1.3 x 10 8 * 
5.66 x 10 6 4.6 x 10 7 4.9 x 10 8 



Factor 0.94 2.08 



2.45 



2.62 



2.65 



2.77 



2.83 



1 The run-times were obtained using Delphi implementations of Ai and A2 on an Intel Core 
Duo CPU E8600 at 3.33 GHz. 

2 Here and in the sequel the algorithm is actually applied to the negative of the intersection 
matrix in question. 
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Clearly, with growing matrix dimension algorithm A 2 shows its advantage over al- 
gorithm A\. (As for the asterisk in the last column: The factors in the fourth line 
of the table seem to suggest a factor of about 3. If this assumption is true, then 
A\ should for n = 64 have a run-time of about 15 x 10 8 ms, which is around 3 
weeks. In our attempts to verify the assumption, we were however not able to get 
results within that amount of time, and we found it technically difficult to check for 
run-times beyond a month.) 

The potential of A2 is demonstrated even more clearly in the case of matrices 
with large definite principal submatrices, e.g., matrices which are negative definite 
themselves. We use a symmetric integer matrix with negative entries —2 on the 
main diagonal, super- and subdiagonal entries 1, and remaining entries 0. Note that 
a similar matrix can in fact be realized as an intersection matrix: For any number 
k consider a surface of the type considered in section [2] of degree k, e.g., the Fermat 
surface of degree k. Among the k 2 lines of the first type we can pick n = 2k A- 1 lines 
with the desired configuration. However, the self-intersections of these lines (and 
therefore the diagonal entries of the matrix) are then 2 — k. 



n 



15 



21 



23 



27 



33 



A 2 



2.73 x 10 2 3.84 x 10 4 4.26 x 10 5 4.5 x 10 6 4.74 x 10 8 
1.1 x 10 2 1.17 x 10 4 1.18 x 10 5 1.23 x 10 6 1.06 x 10 s 



Factor 



2.48 



3.28 



3.61 



3.66 



4.48 
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